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We describe a purely algebraic method for finding the best separable approximation Q to a 
mixed state of a composite 2x2 quantum system, consisting of a decomposition of the state into a 
linear combination of a mixed separable part and a pure entangled one. We prove that, in a generic 
case, the weight of the pure part in the decomposition equals the concurrence of the state. 
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I. INTRODUCTION 



One of the main interests of quantum information theory concerns the nonclassical features connected with the 
nonseparability (or entanglement) of states of composite systems. Since entanglement plays a crucial role in various 
applications in quantum information processing, the problem of characterization of entangled states is of a paramount 
importance. 

It is easy to check whether a pure state of a composite system is separable or entangled. The situation complicates 
in the case of mixed states. A simple and practical necessary criterium of separability is known, but there are no 
known sufficient conditions for higher dimensional composite systems. 

Recently H] , a very interesting description of entanglement was achieved by defining the best separable approxi- 
mation (BSA) of a mixed entangled state. In the simplest case of a 2 x 2 dimensional composite system, it consists 
of a decomposition of the state into a linear combination of a mixed separable part and a pure entangled one. In this 
way, the whole nonseparability properties are concentrated in the pure part. It also provides a natural measure of 
entanglement given by the entanglement of the pure part (well defined for the pure states) multiplied by the weight 
of the pure part in the decomposition. 

In the original paper (T|, the authors proposed a numerical method for finding the BSA in 2 x 2 systems. Some 
analytical results for special states were found in In this paper, we show how to find the BSA of an arbitrary 2x2 
state p in a purely algebraic way, without employing any maximization or optimization procedure. As a byproduct, 
we prove that, in the case that the BSA p s of p is of rank 4, the weight with which the entangled part enters 
the decomposition leading to the BSA, equals another measure of entanglement, namely the concurrence of p ||. 
Furthermore, the pure part is maximally entangled in this case (the last fact was recently proved by other means in 

I). 

The situation is more complicated if the BSA p s is not of full rank. As we will show, for rank(p) = 4 but rank(p s ) < 4 
the components of the BSA are determined by a set of two nonlinear equations which can be easily solved numerically, 
whereas the case of a degenerate p (i.e. rank(p) < 4) can be treated as a limiting case of the full-rank one. It is to stress 
that in these cases there is no simple relation between the concurrence of the state and the weight of the entangled 
part as we were able to prove for rank(p s ) = 4. Presently, we do not have a simple explanation or interpretation of 
this fact which deserves further investigations. 

The paper is organized as follows. In Section |l| we give a maximally shortened account on separability and 



entanglement of mixed states. The main results of the paper are formulated and proved in Section III. The technical 



lemmas used in the proofs of the two main theorems of Section [II are relegated to two Appendices - the first contains 
some more general theorems concerning properties of mixed states of 2 x 2 systems, whereas the second one is mainly 
devoted to a technical lemma concerning relations between spectra of two important matrices obtained from the initial 
mixed state. 

II. SEPARABILITY AND ENTANGLEMENT OF MIXED STATES 

A mixed state p of a bipartite quantum system is separable if it is a convex combination of product states ^ 

k k 
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where pf, pf are legitimate (i.e. hermitian and positive definite) density matrices of the subsystems. 
As observed in 0, a necessary condition for separability of p is that its partial transposition, defined as 



o T ° := $>pf® {pf f, 



(2.2) 



is positive definite, i.e. is als o a l egitimate density matrix for the composite system. (Here, we define the operation of 
partial transposition by Eq (2^2) also in the case of an arbitrary, not necessarily separable state, when pf and pf do 
not need to be positive or/and pi are not all positive - such a decomposition obviously exists for an arbitrary p). Fol- 
low dimensional (2x2 and 2x3) systems the above condition is also sufficient . 

Obviously, the result of partial transposition depends on the basis in subspace TLb . If we change the bases of "Ha 
and TLb by a local transformation U ®V , i.e. by unitary rotations U and V in the spaces TLa and TLb respectively 
(in fact, since the overall phase factor does not play any role, we can assume deiU — 1 = dety, i.e. U, V € SU(2)), 
the matrix p will be transformed according to 



p = U ®Vp{U ®V) ] = ^ Pl Upfu' < ®vpfvt. 



(2.3) 



Consequently, the partial transposition gives 

k 



p 'T B = J2 Pl UpfU^ <g> (Vpf V f ) T = U <g> V* P Tb (U (X) V*)\ 



(2.4) 



where asterix denotes the complex conjugation. From (2.4) it follows that the spectrum of p Ts is basis- independent. 
Observe also the following form of the definition of partial transpose 



(eJ\p\eJ) = {e,r\p T *\e,r), 

where |e, /) denotes the product vector |e) <X> |/). 

In order to quantify the degree of entanglement of two qubit systems, the concurrence was introduced in 

as 

c(p) = max{0, c\ - c 2 — c 3 — c 4 }, 
where c\ > c 2 > c 3 > C4 are the square roots of the (real and positive) eigenvalues of the matrix 

- 







(2.5) 



defined 



(2.6) 



(2.7) 



It is a matter of a straightforward calculation to prove that the concurrence of a pure state, 

\ip) = 01 1 00) +a 2 |Ql) +a 3 |10) +a 4 |ll) = [a u a 2 , a 3 , a 4 ] T , 

equals 

c{ip) = 2\a iai - a 2 a 3 \ = K^|S|V*)I- 



(2.8) 



Due to the normalization condition 1 = (i/>\if>) — \a\\ 2 + |o.2 1 2 + |«3 1 2 + |a.4 1 2 , we have < c(ip) < 1. The maximum 
c(ip) = 1 is attained for the states called maximally entangled. The degree of entanglement (i.e. the concurrence) is 
invariant with respect to local unitary transformations (i.e. transformations of the form U ®V). 

By local transformation, a pure state can be brought to its Schmidt form — Aiei ® fx + \ie 2 ® f 2 , where {ex, e 2 } 
and {/i,/^} are appropriately chosen bases in T-La and TLb- In these bases thus ip — [Ai,0,0, A2] 7, and it is easy to 
show that the most general form of a maximally entangled state in the original bases reads 



|V>) = oi|00) +a 2 |01) T (4110) - fli|H» = 



ax 



"2 



(2.9) 
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III. BEST SEPARABLE APPROXIMATION 



Let p be a generic density matrix for a two qubit system, i.e. a strictly positive definite (i.e. rank 4) 4 x 4 hermitian 
matrix of unit trace. According to jjj, p has a unique decomposition of the form: 

P = (1 - A) \ip) (ip\ + Xp s , (3.1) 

where p s is a separable density ma trix , \tp) is a pure entangled state, and the parameter A € [0, 1] is maximal. In 



the following, we will refer to Eq. as the optimal decomposition of p. The separable part p s is called the best 

separable approximation (BSA) of p, and A its separability. In this section, we will prove the following 
Theorem 1. Let p be an entangled state with rank(p) =4. p = (1 — A) + Xp s is the optimal decomposition 

of p if and only if: 

rank( j o^ B ) = 3, i.e. 3^ p^ B \(f>) — 0, and either 

(i) 3 Q>0 m^\] TB |V> = -c#>, or 

(ii) rank(p s ) = 3, i.e. 3 |<?) p s \4>) = 0, and 3 a ,„> [i#)<0l + [|0)(0|f B ] |V) = 

According to Lemma 2 of Appendix A, |?/>) is maximally entangled in case (i). 

The first condition, rank(p^ B ) = 3 simply states that the BSA p s lies on the boundary between the set of separable 
and the set of entangled states, whereas conditions (i) and (ii) describe the relation between the entangled and 
separable part of the optimal decomposition. Remarkably, the only relevant properties of p s are the vectors \<f>), and 
possibly |</>), in the kernels of p T s B and p s . 

Theorem 1 allows us to check immediately if a given decomposition of p is the optimal one. It also simplifies the 
construction of the BSA for a given p. Indeed, case (i), i.e. any BSA with rank 4, can be solved explicitly, according 
to the following 

Theorem 2. If the BSA p s of p has full rank, the vector \4>) in the (one-dimensional, see Theorem 1) kernel of pJ B 
is an eigenvector of 

Y = Ep TA Ep TB (3.2) 

belonging to the smallest eigenvalue 7 of Y. The weight 1 — A of the entangled part in the optimal decomposition is 
given by 1 — A = 2^/7 = c(p), where c(p) is the concurrence of p. 

Theorem 2 provides a connection between the BSA and the concurrence of p, which was originally || introduced 
as a n auxiliary quantity in order to calculate the entanglement of formation II. Apart from the explicit formula 
( [2.6|) , the concurrence of a mixed state is defined as the minimum of the average concurrence (c) = J2iPi c (' l Pi) over 
all decompositions p = X^PilV'iM^il 01 P m to pure states. After decomposing p s into product states, also the BSA, 
Eq. (|3.lD, defines a particular decomposition, and it follows that 

c(p) < (1 - A)c(V). (3.3) 

This in equa lity implies c(p) + A < 1, which has already been conjectured in p]. According to Theorem 2, e qual ity 



in Eq. (3J3) holds if the BSA of p has full rank (in this case, c{ip) = 1). In other words: the decomposition (3.1) is 



also optimal in the sense that it minimizes the average concurrence. One might assume that t his is true in general 



i.e. also in the second case, rank(,o s ) = 3. Indeed, there exist examples where the inequality (3.3) is saturated also 



in this case, e.g. the generalized Werner states p = x\(j))((f>\ + ^^pl, with \<p) not maximally entangled. (The optimal 



decomposition of these states is given in |2|].) In general, however, we have found that the equality in (3.3) does 
not always hold. Hence, the concurrence of p and the quantity (1 — X)c(ip) provide two inequivalent measures of 
entanglement. Indeed, it is straightforward to show that (1 — A)c(V>) really is a good measure of entanglement, i.e. 
it fulfills the following three conditions ||lo|| : it vanishes if and only if p is separable, is invariant under local unitary 
operations and its expectation value is non- increasing under general local operations (see [ [To[ for details) . 

Before we present the proofs of Theorem 1 and 2, we want to demonstrate how to use the above results in order to 
construct the BSA for a given entangled p of rank 4: first , we calculate the smallest eigenvalue 7 and the corresponding 
eigenvector |</>) of the 4x4 matrix Y, given by Eq. ( |3.2| ). (The eigenvalue 7 is not degenerate, see Lemma 7.) Then, 
we obtain \ip) from condition (i) of Theorem 1, A = 1 — 2^/7, and p s from Eq. ( |3.l[ ). If p s is positive and separable, it 
is the BSA according to Theorem 1. (It is not necessary to check p^ B \ <j>) =0, since this follows from the construction 
of \(f>), see Lemma 7 and Lemmas 3-5.) If not, the BSA has rank 3, and from the second case of Theorem 1 we obtain 
the following set of equations: 



3 



[p\4>)(4>\p] TB 



{iW)p TB \ 

-ap\4>). 



(3.4) 
(3.5) 



Here, we used \ip) — p|0) and (1 — A) -1 = (0|p|0), see Eqs. (3.21) and ( |3.22| ) below. These equations can be solved 
numerically for |0), |0), a, and v. Possibly, ther e exist several solutions, but only one with a, v > and which yields 
a positive and separable state p s via Eq. (3.1). Thereby, we have found the BSA of p in a purely algebraic way, 
without employing any maximization or optimization procedure. 

Finally, we want to show how far this method can be used if rank(p) < 4, in particular if rank(p) = 3, since the 
case rank(p) = 2 has already been solved analytically 0. First, we note that any density matrix p can be obtained 
as a limit from the case of full rank. Thereby, we obtain the following limiting case of Theorem 1 (the complete proof 
will be given below): 

Corollary. Let p be an entangled state with rank(p) < 4. p = (1 — A) \ip) (ip\ + Ap s is the optimal decomposition if 
and only if: 



0, c{4>) > 0, 3 U)Ps 



0, and 3 C 



H0)(0| + [|0)(0|] ifl |V) = -a 



(3.6) 



Although Eq. (3.6) is nearly identical to case (ii) of Theorem 1, it may also arise as a limit from case (i) (with 
v = 0). Unlike in Theorem 1, we must explicitly demand that c(</>) > (in order to exclude solutions which would 
not correspond to the optimal decomposition, see proof of the corollary). 

For the solution of Eq. ( |3.6[) , the following observation is helpful (cf. Lemma 1 of [0): if |0) is not in the kernel of 
p, then \ip) — p\cj) ) and A is given by Eq. ( 3.21 ) (well defined, since is in the range of p). Consequently, Eq. (3.6) 
reduces to Eqs. (3.4 3.5), as in the case of full rank (with the additional constraint c(0) > 0) . Hence, the recipe for 



obtaining the BSA of p with rank(p) = 3 is as follows: first, try to find a solution of Eqs. ( 3.4 ,p.5|) with a, v > and 
c(0) > which yieds a positive and separable p s . If such a solution does not exist, we know that \<f>) must be in the 
kernel of p, which uniquely determines |0) (since we assumed rank(p) = 3). Then, Eq. (3.6) can be solved numerically 
for \ip), |0), a, v and A (after inserting |0) and replacing p T s B by p Ts — (1 — ty[\ip) (i>\] TB ) ■ Furthermore, the following 
fact may be useful: if the kernel of p contains a product vector |e, /), then the solution of Eq. ( |3~6| ) fulfills |0) _L |e, /) 
and |0) _L |e, /*) (see proof of the corollary). 



A. Proof of Theorem 1 



Let p s , as given by Eq. (3.1), be the BSA of p. The maximality condition for A and the uniqueness of the BSA § 
imply [§: 

(a) the state p s + e\ifj)(tp\ is non-separable for e > 0, and 

(b) the state p — (1 — A)| , 0')('0'| is either non-separable or non-positive for each t/j' ^ %l>. 

In order to compact the notation we shall use in the following the notation p, = 1 — A. 
According to the Peres-Horodecki criterion of separability Hf7|, condition (a) implies: 

V e > 3|^> (0 e | P T S B + £ U) l&) < 0. (3.7) 



On the other hand, since p s is separable, the same criterion establishes the positivity of pJ B . Thus, from (3.7) and 
the continuity argument, there is such that 

P T S B \4>) = 0- (3.8) 

Since we assumed rank(p) = 4 and the rank of a projection is one, the rank of p s must be at least three. Then 
rank(pJ B ) = 3, as a consequence of Lemma 1 (cf. Appendix A). 
Now we exploit condition (b). Let's consider 

\p' s :=p-pW)W\, (3.9) 

with 

W) - J5^M, (3.io) 
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where (Aip\Aip) = 1 and (V = 0. (Obviously, any pure state can be written in this form.) 
To the two lowest orders in e we have: 

\p' s = X Ps - ,j (e |V) (AV>| + e |AV) (VI + e 2 |AV)(AVl - e 2 |VXVl) • (3.11) 
In the following, we consider separately two cases of different ranks of p. 
(i) rank(p s ) = 4 

Then, for e small enough, p' s is positive definite for each |AV)- According to the optimality condition (b) above, 
p' s must be non-separable, i.e. there exists such |0') that (0'| p'J B \(f>') < 0. 

Since p^ B has rank 3, 

|0') = |0) + |A0) , (3.12) 



with |A0) — + if e — * 0. Now from ( |3.8| ) we obtain, to the first order in e: 

(A0| ApJ B |A0) - M e<0| [|V) (AVI + |AV) (Vf B 10) < 0. (3.13) 

But p s is, by assumption, separable; consequently p^ B is positive definite 

(A0| \p T s B |A0) > 0, (3.14) 



and ( 3.13 ) implies 

(0| [|V) (AVI + |AV) (V|] Ts 10) >0, (3.15) 
which can be equivalently written as 

Tr{[|V) (AVI + |AV) (V|] Tfl 10) (01} > 0. (3.16) 

For arbitrary operators A and B we have Tt:A Tb B = Tr AB Tb , thus from (|3.16D we obtain 

Tr{[|V) (AVI + |AV) (VII [|0) (0|] Tb } > 0. (3.17) 

This, however, is equivalent to 

(AVI [|0) (0|] Tb |V) + (VI [|0) (0|] Tb |AV) > (3.18) 



for all |AV) -L |V)- Since (3.18) is linear in |AV), changing |AV) into — |AV) reverses the inequality, hence in 
fact it must be that 

(AVI [|0) (0|] Tb |V) = 0. (3.19) 

The above equality must be fulfilled by all |AV) -L |V)- This is possible only if |V) is an eigenvector of 
A = [|0) (0|] Ts - (Note that, although Tg and |0) dep end on the local basis of Hb, the operator A is basis- 
independent, i.e. transforms in the usual way, Eq. ( |2.3| ), under local unitary transformations.) To arrive at the 
first case (i) of Theorem 1, it remains to be shown that the sign of the corresponding eigenvalue a is negative. 
This, however, follows from the limit e — * of Eq. (3.7) 

(0| [|V)(V|] Ts 10) <0, (3.20) 

after using again the identity Tr A Tb B — Tr AB Tb . Furthermore, a cannot be zero - otherwise (ac cording to 
Lemma 2) |0) would be a separable, i.e. a product state: = |e, /), and since pJ B |e, /) = (c.f. (p)^)), we have 
p s \e, f*) = 0, which contradicts the assumption rank(p s ) = 4. From Lemma 2 (c.f. Appendix) we infer that 
|V) is maximally entangled. This provides an alternative proof of the fact proved in [|| that if p and p s are of 
maximal rank then |V) in (3.1) is maximally entangled. 
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(ii) Second case: rank(p s ) < 4 

We assumed that p has rank 4, so rank(p s ) = 3. From Lemma 1 in E|, we know that 



1 - A 



Furthermore, it is easy to check that 



(3.21) 



(3.22) 



fulfills p s 



0. Since rank(p s ) = 3, p' s , Eg. ( 3. 11 ), is positive definite if 



(e|V)(AV|+e|AV)(V|+e 2 |AV)(AV|-e 2 |V)(Vl) \$) < 0. 



(3.23) 



Obviously, this condition is fulfil led if \Aip) _L 
such \Aif>) must fulfill Eq. (|3T9|). 

This is equivalent to 



KVI0) follows from Eq. (|3.2l|) .1 Hence (as in case 1), all 



ivkvi-ivkvimi^if iv) = o, 



(3.24) 



where \ip) is defined such that |V) _L |V) and |V) and |V) span the same two-dimensional subspace as \<p) and 
(We assume that |V) ^ |^>); otherwise, pj, is positive for all |AV), and we get the same result as in the first 
case, which below will turn out to be a special case of the result in the second case.) 

We still have to check the c ase |A V) = \tp). Then, it is always possible to multiply \Aip) by a phase factor such 
that p' s is positive, see Eq. ( 3.23 ) in first order of e. This leads us (as in case 1) to Eq. (|7l|). It follows that 



v (</#>(</#) = (vi mm TB iv> 



(3.25) 



with a nonnegative real parameter v. Otherwise, |AV) could be multiplied by a phase factor such that Eq. 
Q3.23D is fulfilled and Eq. not. (Note that ($\4>) $\ip) ^ 0, since (ip\<t>) ^ follows from Eq. ( |3.21| ), and 

(tp\4>) 7^ from the construction of ip.) 



The two conditions Eqs. (3.25,3.24) are equivalent to the following condition: \ip) is an eigenvector of the 
operator 



A = 



T B 



(3.26) 



To complete the first part of the proof of Theorem 1, we will show now that the corresponding eigenvalue a 
cannot be positive. 

As a consequence of Lemma 2, A has at least three nonnegative eigenvalues. However, there is also at least one 
nonpositive eigenvalue. This follows from the existence of a product vector |e, /) € R(p s ) such that |e, /*) S 
R(p^ B ), as shown in [||, which implies (e,/|A|e,/) = 0. Furthermore, A cannot have more than one zero 
eigenvalue: otherwise, \4>) would have to be a product vector (see Lemma 2), and \4>) would be the corresponding 
partially transposed product vector. Hence, \<j)){4>\ and [\(f>) {4>\] Tb would be identical and proportional to A, and 
|V), as an entangled eigenvector of A, would have to be perpendicular to \<p), i.e. rank(p) = 3, which contradicts 
the assumption rank(p) = 4. 

The above considerations about the spectrum of A are useful for the following reason: let us assume that there 
exists an entangled state p' with a' < which has the property that p(x) = xp + (1 — x)p' is entangled for 
x e [0, 1]. (p' may be a stat e wi th BSA of rank 4, for which we have already shown above that a 1 < 0.) Now, 
the optimal decomposition (3.1) - in particular the eigenvalue a(x) - changes smoothly when varying x from 
to 1 (this follows from the uniqueness of the optimal decomposition). Since, as shown above, A (having one 
nonpositive and three nonnegative eigenvalues) cannot have two zero eigenvalues, a crossing of eigenvalues at 
zero is not possible, and a — a(l) < follows from a' — a(0) < 0. 

It remains to be shown that a state p' with the above properties exists. For this purpose, we consider the Werner 
states p' — y|V')(V'l + ^tt^ with maximally entangled |V')- For these states, it has been shown in [|| that 



G 



the pure state in the optimal decomposition equals \ip') and A' = 3(1 — y)/2. It follows that rank(p' s ) = 4, and 
a' < 0, as shown above (first case). Now, we choose as the eigenvector of [|x)(x|] Tfi with negative eigenvalue 
(such an eigenvalue exists according to Lemma 2), where |x) is an entangled pure state with (x\p Tb \x) < 
(exists, since p is entangled). Using (t/)'\ [\x)(x\] Tb W) = (x\ W) (^'W 3 Ix). [t follows that (x\(p') Tb \x) < 
for large enough y, hence also {x\p( x ) Tb \x) < f° r x € [0) 1]> i- e - Pi x ) is entangled. 

Finally, we will prove the rever se d irection of Theorem 1, i.e. that both cases (i) and (ii) are also sufficient for the 
optimality of the decomposition ( |3.l| ). For this purpose, let us assume that there exists another decomposition with 
larger A. Then, because of the convexity of the set of separable states, such a decomposition with larger A also exists 
in the infinitesimal neighborhood of {A, \ip)}- Hence, for each (infinitesimal small) e > 0, there exists A' = A + AA 
(with AA > and AA — > if e — > 0) and \Aip) _L \ip) such that 

X' P ' S = X Ps + AA|V)M - (1 - A') (e \i/>) (A^| + e |A^> ty\ + e 2 |A^)(A^| - e 2 |V>M) (3.27) 

is separable. Now, let us assume that there exists \<j>) with p^ B \4>) — and either condition (i) or (ii) from Theorem 
(i) is fulfilled. In the following, we will show that both (i) or (ii) lead to a contradiction, since either (4>\p'^ B \<fi) < 
or (4>\p' 8 \4>) <0. 

(i) implies (Aip\ [|0)(0|] Ts IV) = 0, (ip\ [\4>){(t>\] TB < 0, and (A^| U) (</>H Tb |A^>) > 0. (The third inequality 



follows from the spectrum of [|</>) (0|] Ts , see Lemma 2.) Inserting into Eq. ( 3.27 ) immediately yields (4>\p'^ B \4>) < 
0. 

(ii) i mplie s (ip\ [\4>){(j)\] TB \tp) = a - v (ip\4>} {4>\if)) and (Atp\ [\4>){(j)\] TB \ip) = -p{Ai()\4>){(f>\ip). Inserting into Eq. 



(pT27b yields: 

<0|pf B |0) = A\a+(l-\')e 2 (a-(3)-v(4>\p' s \4>), (3.28) 

where (3 = (A-0|A|A-0). Since a < and a < (3 (remember that A = u\(j))(<p\ + [\<f>) (<P\] Tb has three nonnegative 
eigenvalues, i.e. a is the smallest eigenvalue of A), it follows that p' s is either non-positive or non-separable. 



B. Proof of Theorem 2 



Let us assume that Eq. (3.1) is the optimal decomposition of p, with p s of rank 4. According to Theorem 1 (and 
Lemma 2), we know that \tp) is maximally entangled, i.e. c(ip) = 1. Hence, we can use Lemma 3 of Appendix A to 
write 



*Ps ~ P 



pim(w B = p 1b -m( -i-mm 



where \ip) is defined by 



(3.29) 



Consequently, for an arbitrary \<p') 



o < H^W) = (4>'\p Tb W) + Ml W)l 2 - f • 



For \<t>') = |0), the above equation, due to (|3.8|), reads 



0= Mp^^+pM^I 2 -^. 



(3.30) 



(3.31) 



Observe now that because of (i) (Theorem 1) and ( 3.29 ), we can apply Lem ma 5, concluding that \<f>) and \ip) have a 
common Schmidt basis, hence, according to Lemma 4 we can rewrite (3.31) as 



=(01^10) + £<#). 



(3.32) 
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Using the results of the same Lemma we can estimate the last two terms on the right-hand side of (3.3C) by £c(<fc'): 



0<{4>'\ p Tb \^ 



(3.33) 



In o rder to simplify equat ions we are going to consider below, let us make the following observati on. Both equation 
( |3.32 ) and inequality ( 3.33 ) are bilinear in \<j>) if only we calculate the concurrence according to (^8) regardless of 
the normalization of |?/>). Obviously such a quantity is not limited from above, but this will not play any role in the 
following. The final formula will involve only normalized vectors. 

Substituting | (/>') = \ 
in e 



e|A0) (with arbitrary e and |A^>)) to (3.33) and using (3.32), we obtain in the lowest order 



O< £ ((A^|^) + (^|A0) + S dc(0 + eA0) 



de 



e=0 



(3.34) 



From the definition of concurrence ( |2.8| ) we obtain 

c(0 + eA0) = c((f>)+e({A<f>m 

after adjusting the phase of \4>) to make (4>\T,\cf)*) real and positive and using (A0|S|0* 
((f)* |E| Ac/))* which is a consequence of £ = E^ = £*. Thus, we can rewrite (3.34) as 

M 
2 



9 T -|A^) + £((A0|£|0* 



|E|A0)) + o(e 2 



|£|A0» >o, 



|S|A0*) = (A0*|S| 



(3.35) 



valid for an arbitrary |A0). Again, considering ( 3.35 ) for \A<fi) and — |A0), we conclude that in fact ( 3.35 ) is an 
equality 

(A<f>\V) + (*|A0) = 0, 

where 

I*) = p Tb \4>) 

Since |A0) is arbitrary, we have 1^) = and, consequently 



E|</,*). 



Short manipulations using £ = 1 allow for rewriting (3.36) as an eigenvalue equation 

,,2 

E(p TB )*Ep TB | 



4 



(3.36) 



(3.37) 



In Appendix B (Lemma 6), we show that the smallest eigenvalue 7 of Y = E(p Tfi )*Ep TB is given by 7 = c 2 (p)/4, 
where c(p) is the concurrence of p. Furthermore, it follows from Lemma 7 that /i 2 /4 is the smallest eigenvalue of Y, 
since {<p\p TB \4>) < according to Eq. ( |3.36| ). □ 



C. Proof of the Corollary 



Let p be an entangled state with rank(p) < 4, and p = (1 — A) \ip) (ip\ + Xp s its optimal decomposition. Furthermore, 
we define p e := (1 — e)p + el/4 (where I is the 4x4 identity operator). Obviously, rank(p £ ) = 4 and p = lim e ^o p t . 
Since the optimal decomposition of p t varies continously with e, it follows from Theorem 1 in the limit e — > that 



= 0, 



— 0, and 3 C 



H<^| + [k)«>|f B \i>) = -a\i>) 



(3.38) 



This equation includes both cases (i) and (ii) of Theorem 1. (In the former case, v — and \<fi) is an elemen t of the 
kernel of p.) On the other hand, if we assume that |<?!>) is not a product vector (i.e. c(<fi) > 0), then Eq. ( 3.38 ) is 
sufficient for the optimality of the decomposition. This can be proven in the same way as above in the proof of the 
reverse direction of Theorem 1. Indeed, no assumption about the rank of p is needed there, except for showing that 
is not a product vector (which ensures that —a is strictly the smallest eigenvalue of A = v\(f>)((f>\ + [|0)(</>|] Ts ). 
It remains to be shown that c(<p) > is necessary for the optimality of the decomposition. For this purpose, let us 
assume that \4> ) = e, /) is a product vector. From the definition of partial transposition, we know that p s \e, /*) = 0, 
and from Eq. ( J3.38 ) with a > that \tp) _L |e, /*). This implies p\e, /*) = 0. Then, it is easy to show that the BSA of 
p s := (1 - S)p + 5\e,f*)(e,f*\ is given by p Ss = (I - 5)p s + S\e, f*)(e, f*\. Hence, any vector \(f> s ) with ps^ B \(ps) = 
fulfills \(j>g) _L |e, /), which, in the limit S — > 0, contradicts the assumption \<f)) — |e, /). □ 
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APPENDIX A: 

In the Appendix we formulate and prove several Lemmas used in the main part of the paper. 
Lemma 1. Let p s be a two qubit density matrix. If rank( / o^ B ) < 2, then rank(p s ) = rank(pj' f3 ). 
Proof: Since every two-dimensional subspace contains a pro duct vector also the kernel of pj n must do so, i.e. 
p TB \e,f) = 0. It follows that p B \e, /*) = 0. Indeed, from (§]|) we have (e, f*\p\e, /*) = (e, f\p Ts |e, /) = 0, and since 
p s as a density matrix is positive definite, p\e, /*) = 0. By local unitary transformations in both subspaces we can 
choose |e, /*) = |0,0) = |e, /). Equations p\0, 0) = and p Tu |0, 0) = together with hermiticity of both matrices leave 
only six nonvanishing elements in each of them, and by inspection one checks that their characteristic polynomials 
(hence also the spectra) are identical. 

Lemma 2. For an arbitrary \tp} the matrix [|</>)(0|] Tb has eigenvalues 



c c 1 
2' 2'~ 



1 + Vl - c 2 



2 ' 2 

where c = c(0) is the concurrence of |</>). If c > 0, the eigenvector belonging to the negative eigenvalue is maximally 
entangled. 

Proof: The first part of the Lemma is proven by an explicit calculation. In order to prove the second statement, let 
L = U <£> V be a local transformation, and \cf>') = L\4>). Then 



[W){4> 



L' 



where L' — U ®V* . Observe that L 1 is a local transformation, hence it does not influence the concurrence of vectors. 
Now, 



Tfll 



=(0), 



^[W)(4>'\Y B W 



<4>) 



2 UT ' " T " IT ' 2 

where \t/j') = L'\ip). Let us now choose L such that it brings \4>) to its Schmidt basis: 



(Al) 



LI 



Ai 



A 2 



It is now straightforward to show that in (Al) has the form 



V2 







(A2) 



hence \if>') is maximally entangled and the same is true about \t/j) which is obtained from | tp') by a local transformation 
L' . 

(Similar versions of Lemma 1 and Lemma 2 can also be found in |fl2[.) 
Lemma 3. If is maximally entangled then 



where I is the 4x4 identity operator and \tp) is the eigenvector of [|V')('0|] Tb with the negative eigenvalue i.e. 

[|V>)(^|] T ^> = -^). 



(A3) 



(A4) 







According to Lemma 2, \ip) is maximally entangled. 

Proof: Since [\ip) (if>\] Tn is Hermitian, it has, in addition to |-04) 
z = 1,2,3 fulfilling, according to Lemma 2 



1 



i = 1,2,3. 



three other orthogonal eigenvectors l^j), 



(A5) 



Using (A5) and (A4) togethe r wi th the orthonormality of the eigenvectors, (ipi\ipj) = <%, i = 1,2,3,4, one sees that 
the actions of both sides of (A3) give the same results on the complete orthonormal set \ipi), i — 1,2,3,4, which 
establishes the (A3) as a matrix equation. 
Lemma 4. For arbitrary \<p), 



max | {<p\ip) | 



1 1 



(A6) 



The maximum is attained if \tp) and |0) have a 



where the maximum is taken over all maximally entangled | 
common Schmidt basis. 

Proof: By a local unitary transformation (which does not change neither |(</>|-0)| 2 nor the entanglements of \<j>) and 
\tp}) we can bring \<j>) to its Schmidt basis: 



Ai 



A 2 



A; > 0, A 2 + \\ = 1. 



Using the general form (2.9) of a maximally entangled state, we conclude that in the new basis 

|(</#)| 2 = | ffll Ai ± «iA 2 | 2 < |«i| 2 (Ai + A 2 ) 2 = | ai | 2 (1 + 2A X A 2 ) - K| 2 [l + c(0)], 
and the maximum is attained if |ai| 2 is maximal, i.e. |ai| 2 = 1/2 and a 2 = 0, which completes the proof. 



Lemma 5. Let \<p) be an entangled state and \ip) the eigenvector of 

[|0W|] Ts M = -^M, 



Tb with the negative eigenvalue i.e. 



(A7) 



Then \<f>) and \if>) have a common Schmidt basis, where \ip) is the eigenvector of [\ijj) (ip\] TB with the negative eigenvalue, 
i.e. 



i 



Proof: From (A7) we have 

c(4>) 



M U)m TB W) = Tr([|0>(0|] T ^}(^|) = Tr (|0X0| [\^\] 



From Lemma 2 we know that \ip) is maximally entangled. Thus, according to Lemma 3, in the last term we can 
substitute [|V')(V'|] Tfl by \4 ) ){^l>\i consequently: 

hence, from Lemma 4, \<f>) and the maximally entangled \tp) have a common Schmidt basis. 



APPENDIX B: 



Lemma 6. If p is an entangled state, i.e. its concurrence c(p) is positive, then c 2 (p)/4 equals the smallest eigenvalue 
of Y = E(p T *)*Ep T *. 

Proof: If df/4, . . . , d|/4 are the eig enva lues of Y = S( j o Tj3 )*Ep TB and c 2 > . . . > c| the (real and positive, see 
eigenvalues of X = £/9*Ep (c.f. Eq. (2.7)), the following relation holds: 
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d\ = (ci + c 2 + c 3 - c 4 ) 2 , 
d\ = (ci + c 2 - c 3 + c 4 ) 2 , 
dl = (ci - c 2 + c 3 + c 4 ) 2 , 

^4 = ( -C 1 + c 2 + C3 + C 4 ) 2 - 

Indeed, invoking the anticommutation relations for Pauli matrices, we check that for an arbitrary local transforma- 
tion L = U <g) V, U, V £ SU{2) 



If 



(Bl) 



We can thus use local transformations to bring p in X — S / o*E J o and Y = E(/3 Tb )*£ / o Tb to a relatively simple form. 
An arbitrary hermitian p can be decomposed as 



p :-- 



^2 (Ofe^fe ®h + b' k l 2 ® cr fe ) + C„ m er m ® er„, 



(B2) 



with real a' k ,b k , and C mrl . By local transformations, we can bring the 3x3 matrix C to the diagonal form with 
nonnegative diagonal elements p,±,fj,2, and /i 3 ]l3| , |T^ |. The desired transformation changes a' k and b' k to some other 
real a k and b k , hence finally 



P= 4 



a3 + h + P3 bi - ib 2 <xx - ^2 

61 + ib 2 a 3 - 63 - /i 3 /xi + ^ 2 

ai + ia 2 pi + p 2 ~a 3 + b 3 ~ p 3 
Pi — M2 ai + «02 61 + 162 



ai — W2 
-03 - h + M3 



(B3) 



Somewhat tedious but straightforward calculations show that 

Try = TrX 



TrF 4 = TrX 4 -<5 4 



(B4) 
(B5) 
(B6) 
(B7) 



where 



62 
d 2 



6d 
5 



-TrX 2 - -(TrX) 2 , 
2 4 V ; ' 



<5 2 TrX, 



12 
detX 



<S 2 (2TrX 2 + (TrX) 2 



(B8) 

(B9) 

(BIO) 
(Bll) 



On the other hand, as (this time rather short) calculations show, the same relations hold for two diagonal matrices 



X' = diag(c 2 , c 2 ,c 2 ,cl), Y' = diag(d 2 ,d 2 ,d 2 ,dl)/4 



(B12) 



where 



if we choose d = +(detX) 1 / 2 , or 



d 2 
d\ 
dl 
d\ 



d\ = 
d\ = 
d\ = 
dl = 



(cx + c 2 + c 3 - c 4 ) , 
(cx + c 2 - c 3 + c 4 ) 2 , 
(cx - c 2 + c 3 + c 4 ) 2 , 
(-cx + c 2 + c 3 + c 4 ) 2 , 

(-cx + c 2 + c 3 - c 4 ) 2 , 
(-cx + c 2 - c 3 + c 4 ) 2 , 
(-cx - c 2 + c 3 + c 4 ) 2 , 
(cx + c 2 + c 3 + c 4 ) 2 , 



(B13) 
(B14) 
(B15) 
(B16) 



(B17) 
(B18) 
(B19) 
(B20) 
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if d = -{detX) 1 / 2 . Since there is a one to one correspondence between the set of eigenvalues of a n-dimensional 
matrix and the trace s of i ts fir st n powers, t he relation between the eigenvalues cf of X and the eigenvalues df /4 of 
Y must be given by ( B13 - B16 ) or ( B17 - B20 ). The second case, d < 0, is excluded due to the positivity of p. Indeed, 
one checks that: 



6 4 



(TrX) 2 = det(p) > 0. 



(B21) 



The first equality in ( B21 ) follows from (B8) whereas the second is established by an explicit calculation using (B3) 
and the definition of X in terms of p. 



d\ (B13-B16) of AY are real and positive, and the smallest eigenvalue d\ equals c(pY 



Thus the eigenvalues d\ , 
see Eq. (gg). 

Lemma 7. If rank( ( o) > 3, where p is an entangled state, the smallest eigenvalue of Y 



E(p Tfl )*Ep J ' B is non 



degenerate. If |</>4) 
following holds: 



denotes the corresponding eigenvector, and \(/>i), i — 1,2,3, the other three eigenvectors, the 



h\p TB \<j>i) = --c(p)c(^ 4 ), 



\P 



k) > 0, 



1,2,3. 



(B22) 
(B23) 



Proof: Let df/4 > . .. > d|/4 denote the (real and positive, see Lemma 6) eigenvalues of Y, and c\ > ._ L , > cj the 
eigenvalues of X — E/9*Ep. According to Lemma 6, the relation between di and a is given by Eqs. (pnjpH), m 
particular d^ = c(p). From c(p) > and the definition of concurrence, Eq. ( |2.6| ), it follows that C\ > c 2 . Now, i f 
rank(p) > 3, it is easy to show that C2 > (since rank(E) = 4 and therefore rank(X) > 2), and then Eqs. (B15 |B1^ ) 
imply c?4 < d$. Hence, d\/A is a non-degenerate eigenvalue. 

By splitting the eigenvalue equation Y\(f>i) = jdf\(j)i) (with real di) into its real and imaginary part, one can derive 
that \<f>i) fulfills Sp TB |0i 
P Tb , we conclude that 



2 ^ 



), where e lXi is a phase factor. Using E 2 = 1, Eq. (2.8), and the hermiticity of 



P TB \4>i 



(B24) 



In order to complete the proof of Lemma 7, it remains to be shown that the sign on the right hand side must be negative 
for i = 4 and nonnegative for z = 1,2, 3. Because of continuity, it is sufficient to consider the case rank(,o) = 4. Then, 
\(pi) cannot be a produ ct ve ctor (since inserting 1^) = |e,/) into Eq. (B24) would imply (e, f*\p\e, /*) = 0), i.e. the 
right hand side of Eq. (B24) cannot be zero (di > follows from d^ — c (p) > 0). Now, if p is infinitesimally close to an 
entangled pure state, p it is easy to check that, indeed, Eq. ( B24 ) is valid with the minus sign for i = 4 and 

the plus sign for i = 1,2,3. (For \f) = [Ai, 0, 0, A 2 ] T , one finds that |0i, 2 ) = [A 2 , 0, 0, ±Ai] T , \<f> 3 ) = [0, 1, 1, 0] T /V2, 
and 1^4) = [0, 1, — 1, 0] T /\/2-) Next, we consider the one parameter family p(X') = p'\tp)(tp\ + \'p s , with pf = 1 — A' 
and A' € [0, A], where p s is the BSA of p = p(X). Since pJ B is positive, (x\p Tb \x) < implies (x|p(A / ) Ts \x) < 0, hence 
c(p(X')) > for all A' S [0, A]. Finally, continuity implies that the sign of the right hand side of Eq. ( |B24 ) does not 
change when increasing A' from to A. 



M. Lewenstein and A. Sanpera, Phys. Rev. Lett. 80, 2261 (1998). 
B.-G. Englert and N. Metwally, J. Mod. Opt. 47, 2221 (2000). 
W. K. Wootters, Phys. Rev. L ett. 80, 2245 (1998 ). 



S. Karnas and M. Lewenstein, |quant-ph/0011066| v2 (2000) 
R. F. Werner, Phys. Rev. A 40, 4277 (1989). 

A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. A 223, 1 (1996). 
S. K. Hill and W. K. Wootters, Phys. Rev. Lett. 78, 5022 (1997). 

C. H. Bennet, D. P. DiVincenzo, J. Smolin, and W. K. Wootters, Phys. Rev. A 54, 3824 (1996). 
V. Vedral and M. B. Plenio, Phys. Rev. A 57, 1619 (1998). 

B. -G. Englert and N. Metwally, Appl. Phys. B - Lasers O 72, 35 (2000). 
A. Sanpera, R. Tarrach, and G. Vidal, Phys. Rev. A 58, 826 (1998). 



12 



[13] R. Horodecki and M. Horodecki, Phys. Rev. A 54, 1838 (1996). 
[14] M. Kus and K. Zyczkowski, Phys. Rev. A 63, 032307 (2001). 



13 



